Multi-cavity molecular descriptor interconnections: Enhanced protocol for prediction of serum albumin drug binding.

The role of human serum albumin (HSA) in the transport of molecules predicates its involvement in the determination of drug distribution and metabolism. Optimization of ADME properties are analogous to HSA binding thus this is imperative to the drug discovery process. Currently, various in silico predictive tools exist to complement the drug discovery process, however, the prediction of possible ligand-binding sites on HSA has posed several challenges. Herein, we present a strong and deeper-than-surface case for the prediction of HSA-ligand binding sites using multi-cavity molecular descriptors by exploiting all experimentally available and crystallized HSA-bound drugs. Unlike previously proposed models found in literature, we established an in-depth correlation between the physicochemical properties of available crystallized HSA-bound drugs and different HSA binding site characteristics to precisely predict the binding sites of investigational molecules. Molecular descriptors such as the number of hydrogen bond donors (nHD), number of heteroatoms (nHet), topological polar surface area (TPSA), molecular weight (MW), and distribution coefficient (LogD) were correlated against HSA binding site characteristics, including hydrophobicity, hydrophilicity, enclosure, exposure, contact, site volume, and donor/acceptor ratio. Molecular descriptors nHD, TPSA, LogD, nHet, and MW were found to possess the most inherent capacities providing baseline information for the prediction of serum albumin binding site. We believe that these associations may form the bedrock for establishing a solid correlation between the physicochemical properties and Albumin binding site architecture. Information presented in this report would serve as critical in provisions of rational drug designing as well as drug delivery, bioavailability, and pharmacokinetics.

Exploring the potential of biologically active phenolic acids from marine natural products as anticancer agents targeting the epidermal growth factor receptor.

The epidermal growth factor receptor (EGFR) dimerizes upon ligand bindings to the extracellular domain that initiates the downstream signaling cascades and activates intracellular kinase domain. Thus, activation of autophosphorylation through kinase domain results in metastasis, cell proliferation, and angiogenesis. The main objective of this research is to discover more promising anti-cancer lead compound against EGRF from the phenolic acids of marine natural products using in-silico approaches. Phenolic compounds reported from marine sources are reviewed from previous literatures. Furthermore, molecular docking was carried out using the online tool CB-Dock. The molecules with good docking and binding energies scores were subjected to ADME, toxicity and drug-likeness analysis. Subsequently, molecules from the docking experiments were also evaluated using the acute toxicity and MD simulation studies. Fourteen phenolic compounds from the reported literatures were reviewed based on the findings, isolation, characterized and applications. Molecular docking studies proved that the phenolic acids have good binding fitting by forming hydrogen bonds with amino acid residues at the binding site of EGFR. Chlorogenic acid, Chicoric acid and Rosmarinic acid showed the best binding energies score and forming hydrogen bonds with amino acid residues compare to the reference drug Erlotinib. Among these compounds, Rosmarinic acid showed the good pharmacokinetics profiles as well as acute toxicity profile. The MD simulation study further revealed that the lead complex is stable and could be future drug to treat the cancer disease. Furthermore, in a wet lab environment, both in-vitro and in-vivo testing will be employed to validate the existing computational results.Communicated by Ramaswamy H. Sarma.

Selected phytochemicals of Momordica charantia L. as potential anti-DENV-2 through the docking, DFT and molecular dynamic simulation.

Dengue fever is now one of the major global health concerns particularly for tropical and sub-tropical countries. However, there has been no FDA approved medication to treat dengue fever. Researchers are looking into DENV NS5 RdRp protease as a potential therapeutic target for discovering effective anti-dengue agents. The aim of this study to discover dengue virus inhibitor from a set of five compounds from Momordica charantia L. using a series of in-silico approaches. The compounds were docked into the active area of the DENV-2 NS5 RdRp protease to obtain the hit compounds. The successful compounds underwent additional testing for a study on drug-likeness similarity. Our study obtained Momordicoside-I as a lead compound which was further exposed to the Cytochrome P450 (CYP450) toxicity analysis to determine the toxicity based on docking scores and drug-likeness studies. Moreover, DFT studies were carried out to calculate the thermodynamic, molecular orbital and electrostatic potential properties for the lead compound. Moreover, the lead compound was next subjected to molecular dynamic simulation for 200 ns in order to confirm the stability of the docked complex and the binding posture discovered during docking experiment. Overall, the lead compound has demonstrated good medication like qualities, non-toxicity, and significant binding affinity towards the DENV-2 RdRp enzyme.Communicated by Ramaswamy H. Sarma.

Investigation of Simultaneous and Sequential Cooperative Homotropic Inhibitor Binding to the Catalytic Chamber of SARS-CoV-2 RNA-dependent RNA Polymerase (RdRp).

In our previous report, the unique architecture of the catalytic chamber of the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp), which harbours two distinctive binding sites, was fully characterized at molecular level. The significant differences in the two binding sites BS1 and BS2 in terms of binding pockets motif, as well as the preferential affinities of eight anti-viral drugs to each of the two binding sites were described. Recent Cryogenic Electron Microscopy (Cryo-EM) studies on the RdRp revealed that two suramin molecules, a SARS-CoV-2 inhibitor, bind to RdRp in two different sites with distinctive interaction landscape. Here, we provide the first account of investigating the combined inhibitor binding to both binding sites, and whether the binding of two inhibitors molecules concurrently is "Cooperative binding" or not. It should be noted that the binding of inhibitors to different sites do not necessary constitute mutually independent events, therefore, we investigated two scenarios to better understand cooperativity: simultaneous binding and sequential binding. It has been demonstrated by binding free energy calculations (MM/PBSA) and piecewise linear potential (PLP) interaction energy analysis that the co-binding of two suramin molecules is not cooperative in nature; rather, when compared to individual binding, both molecules adversely affect one another's binding affinities. This observation appeared to be primarily due to RdRp's rigidity, which prevented both ligands from fitting comfortably within the catalytic chamber. Instead, the suramin molecules showed a tendency to change their orientation within the binding pockets in order to maintain their binding to the protein, but at the expense of the ligand internal energies. Although co-binding resulted in the loss of several important key interactions, a few interactions were conserved, and these appear to be crucial in preserving the binding of ligands in the active site. The structural and mechanistic details of this study will be useful for future research on creating and developing RdRp inhibitors against SARS-CoV-2.

Therapeutic Path to Triple Knockout: Investigating the Pan-inhibitory Mechanisms of AKT, CDK9, and TNKS2 by a novel 2-Phenylquinazolinone derivative in Cancer Therapy- An In-Silico Investigation.

BACKGROUND: Blocking the oncogenic Wnt//ß-catenin pathway has of late been investigated as a viable therapeutic approach in the treatment of cancer. This involves the multi-targeting of certain members of the tankyrase-kinase family; tankyrase 2 (TNKS2), protein kinase B (AKT), and cyclin-dependent kinase 9 (CDK9), which propagate the oncogenic Wnt/ß-catenin signalling pathway. METHODS: During a recent investigation, the pharmacological activity of 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one was repurposed to serve as a 'triple-target' inhibitor of TNKS2, AKT and CDK9. Yet, the molecular mechanism that surrounds its multi-targeting activity remains unanswered. As such, this study aims to explore the pan-inhibitory mechanism of 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one towards AKT, CDK9, and TNKS2, using in silico techniques. RESULTS: Results revealed favourable binding affinities of -34.17 kcal/mol, -28.74 kcal/mol, and -27.30 kcal/mol for 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one towards TNKS2, CDK9, and AKT, respectively. Pan-inhibitory binding of 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one is illustrated by close interaction with specific residues on tankyrase-kinase. Structurally, 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one had an impact on the flexibility, solvent-accessible surface area, and stability of all three proteins, which was illustrated by numerous modifications observed in the unbound as well as the bound states of the structures, which evidenced the disruption of their biological function. Prediction of the pharmacokinetics and physicochemical properties of 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one further established its inhibitory potential, evidenced by the favourable absorption, metabolism, excretion, and minimal toxicity properties. CONCLUSION: The following structural insights provide a starting point for understanding the pan-inhibitory activity of 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one. Determining the criticality of the interactions that exist between the pyrimidine ring and catalytic residues could offer insight into the structure-based design of innovative tankyrase-kinase inhibitors with enhanced therapeutic effects.

A novel small molecule inhibitor of p38⍺ MAP kinase augments cardiomyocyte cell cycle entry in response to direct cell cycle stimulation.

BACKGROUND AND PURPOSE: Myocardial infarction (MI) is the leading cause of mortality globally due in part to the limited ability of cardiomyocytes (CMs) to regenerate. Recently, we demonstrated that overexpression of four-cell cycle factors, CDK1, CDK4, cyclin B1 and cyclin D1 (4F), induced cell division in ~20% of the post-mitotic CMs overexpressed 4F. The current study aims to identify a small molecule that augments 4F-induced CM cycle induction. EXPERIMENTAL APPROACH, KEY RESULTS: Screening of small molecules with a potential to augment 4F-induced cell-cycle induction in 60-day-old mature human induced pluripotent cardiomyocytes (hiPS-CMs) revealed N-(4,6-Dimethylpyridin-2-yl)-4-(pyridine-4-yl)piperazine-1-carbothioamide (NDPPC), which activates cell cycle progression in 4F-transduced hiPS-CMs. Autodock tool and Autodock vina computational methods showed that NDPPC has a potential interaction with the binding site at the human p38⍺ mitogen-activated protein kinase (p38⍺ MAP kinase), a critical negative regulator of the mammalian cell cycle. A p38 MAP kinase activity assay showed that NDPPC inhibits p38⍺ with 5-10 times lower IC50 compared to the other P38 isoforms in a dose-dependent manner. Overexpression of p38⍺ MAP kinase in CMs inhibited 4F cell cycle induction, and treatment with NDPPC reversed the cell cycle inhibitory effect. CONCLUSION AND IMPLICATIONS: NDPPC is a novel inhibitor for p38 MAP kinase and is a promising drug to augment CM cell cycle response to the 4F. NDPPC could become an adjunct treatment with other cell cycle activators for heart failure treatment.

Artificial Intelligence, Machine Learning, and Big Data for Ebola Virus Drug Discovery.

The effect of Ebola virus disease (EVD) is fatal and devastating, necessitating several efforts to identify potent biotherapeutic molecules. This review seeks to provide perspectives on complementing existing work on Ebola virus (EBOV) by discussing the role of machine learning (ML) techniques in the prediction of small molecule inhibitors of EBOV. Different ML algorithms have been used to predict anti-EBOV compounds, including Bayesian, support vector machine, and random forest algorithms, which present strong models with credible outcomes. The use of deep learning models for predicting anti-EBOV molecules is underutilized; therefore, we discuss how such models could be leveraged to develop fast, efficient, robust, and novel algorithms to aid in the discovery of anti-EBOV drugs. We further discuss the deep neural network as a plausible ML algorithm for predicting anti-EBOV compounds. We also summarize the plethora of data sources necessary for ML predictions in the form of systematic and comprehensive high-dimensional data. With ongoing efforts to eradicate EVD, the application of artificial intelligence-based ML to EBOV drug discovery research can promote data-driven decision making and may help to reduce the high attrition rates of compounds in the drug development pipeline.

In silico evaluation of usnic acid derivatives to discover potential antibacterial drugs against DNA gyrase B and DNA topoisomerase IV.

Due to the rising increase in infectious diseases brought on by bacteria and anti-bacterial drug resistance, antibacterial therapy has become difficult. The majority of first-line antibiotics are no longer effective against numerous germs, posing a new hazard to global human health in the 21st century. Through the drug-likeness screening, 184 usnic acid derivatives were selected from an in-house database of 340 usnic acid compounds. The pharmacokinetics (ADMET) prediction produced fifteen hit compounds, of which the lead molecule was subsequently obtained through a molecular docking investigation. The lead compounds, labelled compound-277 and compound-276, respectively, with the substantial binding affinity towards the enzymes were obtained through further docking simulation on the DNA gyrase and DNA topoisomerase proteins. Additionally, molecular dynamic (MD) simulation was performed for 300 ns on the lead compounds in order to confirm the stability of the docked complexes and the binding pose discovered during docking tests. Due to their intriguing pharmacological characteristics, these substances may be promising therapeutic candidate for anti-bacterial medication.Communicated by Ramaswamy H. Sarma.

Multi-dimensional structural footprint identification for the design of potential scaffolds targeting METTL3 in cancer treatment from natural compounds.

CONTEXT: [Formula: see text]-adenosine-methyltransferase (METTL3) is the catalytic domain of the 'writer' proteins which is involved in the post modifications of [Formula: see text]-methyladinosine ([Formula: see text]). Though its activities are essential in many biological processes, it has been implicated in several types of cancer. Thus, drug developers and researchers are relentlessly in search of small molecule inhibitors that can ameliorate the oncogenic activities of METTL3. Currently, STM2457 is a potent, highly selective inhibitor of METTL3 but is yet to be approved. METHODS: In this study, we employed structure-based virtual screening through consensus docking by using AutoDock Vina in PyRx interface and Glide virtual screening workflow of Schrodinger Glide. Thermodynamics via MM-PBSA calculations was further used to rank the compounds based on their total free binding energies. All atom molecular dynamics simulations were performed using AMBER 18 package. FF14SB force fields and Antechamber were used to parameterize the protein and compounds respectively. Post analysis of generated trajectories was analyzed with CPPTRAJ and PTRAJ modules incorporated in the AMBER package while Discovery studio and UCSF Chimera were used for visualization, and origin data tool used to plot all graphs. RESULTS: Three compounds with total free binding energies higher than STM2457 were selected for extended molecular dynamics simulations. The compounds, SANCDB0370, SANCDB0867, and SANCDB1033, exhibited stability and deeper penetration into the hydrophobic core of the protein. They engaged in relatively stronger intermolecular interactions involving hydrogen bonds with resultant increase in stability, reduced flexibility, and decrease in the surface area of the protein available for solvent interactions suggesting an induced folding of the catalytic domain. Furthermore, in silico pharmacokinetics and physicochemical analysis of the compounds revealed good properties suggesting these compounds could serve as promising MEETL3 entry inhibitors upon modifications and optimizations as presented by natural compounds. Further biochemical testing and experimentations would aid in the discovery of effective inhibitors against the berserk activities of METTL3.

The Unusual Architecture of RNA-Dependent RNA Polymerase (RdRp)'s Catalytic Chamber Provides a Potential Strategy for Combination Therapy against COVID-19.

The unusual and interesting architecture of the catalytic chamber of the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) was recently explored using Cryogenic Electron Microscopy (Cryo-EM), which revealed the presence of two distinctive binding cavities within the catalytic chamber. In this report, first, we mapped out and fully characterized the variations between the two binding sites, BS1 and BS2, for significant differences in their amino acid architecture, size, volume, and hydrophobicity. This was followed by investigating the preferential binding of eight antiviral agents to each of the two binding sites, BS1 and BS2, to understand the fundamental factors that govern the preferential binding of each drug to each binding site. Results showed that, in general, hydrophobic drugs, such as remdesivir and sofosbuvir, bind better to both binding sites than relatively less hydrophobic drugs, such as alovudine, molnupiravir, zidovudine, favilavir, and ribavirin. However, suramin, which is a highly hydrophobic drug, unexpectedly showed overall weaker binding affinities in both binding sites when compared to other drugs. This unexpected observation may be attributed to its high binding solvation energy, which disfavors overall binding of suramin in both binding sites. On the other hand, hydrophobic drugs displayed higher binding affinities towards BS1 due to its higher hydrophobic architecture when compared to BS2, while less hydrophobic drugs did not show a significant difference in binding affinities in both binding sites. Analysis of binding energy contributions revealed that the most favorable components are the ΔEele, ΔEvdw, and ΔGgas, whereas ΔGsol was unfavorable. The ΔEele and ΔGgas for hydrophobic drugs were enough to balance the unfavorable ΔGsol, leaving the ΔEvdw to be the most determining factor of the total binding energy. The information presented in this report will provide guidelines for tailoring SARS-CoV-2 inhibitors with enhanced binding profiles.

Intermolecular And Dynamic Investigation of The Mechanism of Action of Reldesemtiv on Fast Skeletal Muscle Troponin Complex Toward the Treatment of Impaired Muscle Function.

Muscle weakness as a secondary feature of attenuated neuronal input often leads to disability and sometimes death in patients with neurogenic neuromuscular diseases. These impaired muscle function has been observed in several diseases including amyotrophic lateral sclerosis, Charcot-Marie-Tooth, spinal muscular atrophy and Myasthenia gravis. This has spurred the search for small molecules which could activate fast skeletal muscle troponin complex as a means to increase muscle strength. Discovered small molecules have however been punctuated by off-target and side effects leading to the development of the second-generation small molecule, Reldesemtiv. In this study, we investigated the impact of Reldesemtiv binding to the fast skeletal troponin complex and the molecular determinants that condition the therapeutic prowess of Redesemtiv through computational techniques. It was revealed that Reldesemtiv binding possibly potentiates troponin C compacting characterized by reduced exposure to solvent molecules which could favor the slow release of calcium ions and the resultant sensitization of the subunit to calcium. These conformational changes were underscored by conventional and carbon hydrogen bonds, pi-alkyl, pi-sulfur and halogen interactions between Reldesemtiv the binding site residues. Arg113 (-3.96 kcal/mol), Met116 (-2.23 kcal/mol), Val114 (-1.28 kcal/mol) and Met121 (-0.63 kcal/mol) of the switch region of the inhibitory subunit were among the residues that contributed the most to the total free binding energy of Reldesemtiv highlighting their importance. These findings present useful insights which could lay the foundation for the development of fast skeletal muscle small molecule activators with high specificity and potency.

Chetomin, a SARS-CoV-2 3C-like Protease (3CLpro) Inhibitor: In Silico Screening, Enzyme Docking, Molecular Dynamics and Pharmacokinetics Analysis.

The emergence of the Coronavirus Disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has led to over 6 million deaths. The 3C-like protease (3CLpro) enzyme of the SARS-CoV-2 virus is an attractive druggable target for exploring therapeutic drug candidates to combat COVID-19 due to its key function in viral replication. Marine natural products (MNPs) have attracted considerable attention as alternative sources of antiviral drug candidates. In looking for potential 3CLpro inhibitors, the MNP database (>14,000 molecules) was virtually screened against 3CLpro with the assistance of molecular docking computations. The performance of AutoDock and OEDocking software in anticipating the ligand-3CLpro binding mode was first validated according to the available experimental data. Based on the docking scores, the most potent MNPs were further subjected to molecular dynamics (MD) simulations, and the binding affinities of those molecules were computed using the MM-GBSA approach. According to MM-GBSA//200 ns MD simulations, chetomin (UMHMNP1403367) exhibited a higher binding affinity against 3CLpro than XF7, with ΔGbinding values of −55.5 and −43.7 kcal/mol, respectively. The steadiness and tightness of chetomin with 3CLpro were evaluated, revealing the high stabilization of chetomin (UMHMNP1403367) inside the binding pocket of 3CLpro throughout 200 ns MD simulations. The physicochemical and pharmacokinetic features of chetomin were also predicted, and the oral bioavailability of chetomin was demonstrated. Furthermore, the potentiality of chetomin analogues −namely, chetomin A-D− as 3CLpro inhibitors was investigated. These results warrant further in vivo and in vitro assays of chetomin (UMHMNP1403367) as a promising anti-COVID-19 drug candidate.

Probing into the Flap-dimer Dynamics of the Mycobacterium tuberculosis Kasa Enzyme Binding Landscape Provides the Underlying Inhibitory Mechanisms of JSF-3285 and 5G.

BACKGROUND: ß-ketoacyl-ACP synthase I (KasA I) enzyme is crucial in mycolic acid synthesis via catalytic condensation reactions, hence implicated in M. tuberculosis's virulence and drug resistance. Presently, there is no known potent KasA inhibitor; thiolactomycin lacks potency. Recently reported indazole compounds JSF-3285/tr1DG167 and 5G/tr2DG167 inhibit the KasA through binding to the substrate cavity. However, the molecular mechanism is still unclear, and the unknown resistance mechanisms raise concerns about JSF-3285's novelty. METHODS: This study is the first to report the flap dimer opening and closing of the KasA pocket using combined metrics to define the symmetry impact of the flap-dimer motions and investigate the underlying inhibitory mechanism of tr1DG167 andtr2DG167 using all-atom MD simulation. RESULTS: The distance/d1 between the flap (PRO147) and dimer (LEU205) residues; TriC-α angle (θ1: PRO147-VAL83-LEU205 & θ2: PRO147-GLU199-LEU205); and the dihedral angle (Φ) were applied to investigate the flap "twisting" and dimer shift closing due to concerted motion by adjacent glycine-rich and glutamic acid-rich loops around the active site during the binding pocket's opening. The full flap-dimer of the unbound opens at 230 ns (d1 = 21.51 Å), corresponding to the largest TriC-α angle θ1 44.5° as θ2 is unreliable to describe the flap-dimer motion. The overall averages θ1 and θ2 for the bounds were ~23.13° and ~23.31°, respectively. Thus, the degree of KasA flap dimer opening is best investigated by distance and θ1. BFE (Kcal/mol) of -44.05 (tr1DG167) showed a higher affinity for the pocket than tr2DG167-KasA (-32.16). Both tr1DG167 and tr2DG167 formed hydrophobic interactions with LEU116, GLY117, ALA119, and tr1DG167 formed strong H-bonds with GLU199. The average RMSD of 2.80 Å (Apo) and RoG of 20.97 Å showed that KasA is less stable and less tightly packed without the inhibitors. CONCLUSION: These findings provide a background for a new structure-based design of novel KasA inhibitors.

Pharmacophore-based virtual screening and in-silico study of natural products as potential DENV-2 RdRp inhibitors.

Dengue fever is a significant public health concern throughout the world, causing an estimated 500,000 hospitalizations and 20,000 deaths each year, despite the lack of effective therapies. The DENV-2 RdRp has been identified as a potential target for the development of new and effective dengue therapies. This research's primary objective was to discover an anti-DENV inhibitor using in silico ligand- and structure-based approaches. To begin, a ligand-based pharmacophore model was developed, and 130 distinct natural products (NPs) were screened. Docking of the pharmacophore-matched compounds were performed to the active site of DENV-2 RdRp protease . Eleven compounds were identified as potential DENV-2 RdRp inhibitors based on docking energy and binding interactions. ADMET and drug-likeness were done to predict their pharmacologic, pharmacokinetic, and drug-likeproperties . Compounds ranked highest in terms of pharmacokinetics and drug-like appearances were then subjected to additional toxicity testing to determine the leading compound. Additionally, MD simulation of the lead compound was performed to confirm the docked complex's stability and the binding site determined by docking. As a result, the lead compound (compound-108) demonstrated an excellent match to the pharmacophore, a strong binding contact and affinity for the RdRp enzyme, favourable pharmacokinetics, and drug-like characteristics. In summary, the lead compound identified in this study could be a possible DENV-2 RdRp inhibitor that may be further studied on in vitro and in vivo models to develop as a drug candidate.Communicated by Ramaswamy H. Sarma.

Investigating the Impact of Covalent and Non-covalent Binding Modes of Inhibitors on Bruton's Tyrosine Kinase in the Treatment of B Cell Malignancies - Computational Insights.

BACKGROUND: Bruton tyrosine kinase plays a key role in the survival, proliferation, activation, and differentiation of B-lineage cells and the signaling of other receptors. It is overexpressed and constitutively active in the pathogenesis of B cell malignancies and has therefore become a target for therapeutic intervention. Some success has been achieved in the discovery of small molecules, especially in the development of irreversible inhibitors. However, these inhibitors are punctuated by off target effects and have also become less effective in patients with mutations at Cys481. This motivated the search for inhibitors with improved efficacy and different binding modes. METHODS: In this study, we employed two new second generation inhibitors with different binding modes, Zanubrutinib and AS-1763, which are at various levels of clinical trials, to highlight the molecular determinants in the therapeutic inhibition of BTK through computational studies. RESULTS: This study revealed that Zanubrutinib and AS-1763 exhibited free total binding energies of -98.76 ± 4.63 kcal/mol and -51.81 ± 9.94 kcal/mol, respectively, with Zanubrutinib engaging in peculiar hydrogen bond interactions with the hinge residues Glu475 and Met477 including Asn484 and Tyr485 while AS-1763 engaged Lys430, Asp539, and Arg525. These residues contributed the most towards the free total binding energy with energies above -1.0 kcal/mol. The compounds further interacted differentially with other binding site residues through pi-alkyl, pi-cation, pianion, pi-pi-T-shaped, pi-sigma, pi-sulfur and pi-donor hydrogen bonds, and Van der Waals interactions. These interactions resulted in differential fluctuations of the residues with the consequential unfolding of the protein. CONCLUSION: Insights herein would be useful in guiding the discovery of more selective and potent small molecules.

Deciphering the molecular mechanisms of selective non-covalency demonstrated differentially by 9-Allylnaphtho[1,8-ef]isoindole-7,8,10(9H)-trione (C11) against fibroblast growth factor receptors 1-4.

The specific inhibition of aberrant Fibroblast Growth Factor Receptors (FGFRs) has been identified as a feasible strategy to therapeutically ameliorate their respective carcinogenic involvements. High homology among these proteins has however limited efforts towards the discovery of selective small-molecule compounds due to undesirable effects elicited by pan-FGFR inhibitors. A recent study showed the selective activity of a new compound C11 which was >52 times more potent against FGFR1 than FGFR2 and FGFR3, and 4 times than FGFR4. This C11 selective non-covalency was investigated in this study using computational methods since it has remained unresolved. Structural findings revealed that C11 enhanced structural perturbations in FGFR1 with less prominent effects in other FGFRs. High deviations also characterized the C11-bound active pocket of FGFR1 with notable fluctuations across the constituent P-loop, αC helix, hinge region, catalytic, and activation loops. These induced motions were essential for optimal C11 motion an d positioning of its phenalenone ring and prop-2-en-l-yl moiety at the FGFR1 active pocket to interact stably and strongly with A564FGFR1, L484FGFR1, Y563FGFR1, and E562FGFR1 which as well had high energy contributions. C11 exhibited highly unstable binding in F GFRs2-3 with a more steady interaction with FGFR4. Free binding energy (ΔGbind) analyses further estimated the highest interaction energy for C11-FGFR1 with favorable desolvation energy that indicated a deep hydrophobic pocket binding for C11 in FGFR1 compared to other FGFRs. We believe rational insights from this study will contribute to the structure-based design of highly specific FGFR1 inhibitors.Communicated by Ramaswamy H. Sarma.

Impact of compound mutations I1171N + F1174I and I1171N + L1198H on the structure of ALK in NSCLC pathogenesis: atomistic insights.

Anaplastic lymphoma kinase (ALK) fusion genes are found in 3%-5% of non-small cell lung cancers (NSCLCs). NSCLC is the most common type of lung cancer, accounting for 84% of all lung cancer diagnoses. Available treatment options for ALK-positive NSCLCs involve the use of ALK tyrosine kinase inhibitors (ALK-TKIs) which have shown to be effective with a high response rate. Nonetheless, the emergence of multiple compound mutations such as I1171N + F1174I or I1171N + L1198H has been reported to cause resistance to all approved ALK-TKIs. However, the underlying molecular mechanisms surrounding the impact of these compound mutants remain poorly understood. Hence, we performed molecular dynamics simulations to characterize the structural effects and functional implications of these compound mutations. Findings revealed a destabilizing effect on ALK by mutants as compared to the wild-type ALK structure. Also, further insights revealed a lower root-mean-squared fluctuation, radius of gyration, and solvent-accessible surface area values of I1171N + F1174I and I1171N + L1198H ALK compound mutations suggesting that the mutants have a more compact structure and a smaller surface area than the wild-type protein. The mutants also distorted the activation loop residues (Tyr1278, Tyr1282, and Tyr1283) in the ALK structure, which further identify them as possible disruptors of phosphorylation. In contrast to wild conformation, the mutant conformations exhibited a reduced node degree in their residue interaction networks. Collectively, our findings provide deeper insights into the deleterious effects of I1171N + F1174I and I1171N + L1198H ALK compound mutations, which may contribute to NSCLC pathogenesis.Communicated by Ramaswamy H. Sarma.

Inside the cracked kernel: establishing the molecular basis of AMG510 and MRTX849 in destabilising KRASG12C mutant switch I and II in cancer treatment.

The Kirsten rat sarcoma oncoprotein (KRAS) has been punctuated by drug development failures for decades due to frequent mutations that occur mostly at codon 12 and the seemingly intractable targeting of the protein. However, with advances in covalent targeting, the oncoprotein is being expunged from the 'undruggable' list of proteins. This feat has seen some covalent drugs at different stages of clinical trials. The advancement of AMG510 and MRTX849 as inhibitors of cysteine mutated KRAS (KRASG12C) to phase-III clinical trials informed the biased selection of AMG510 and MRTX849 for this study. Despite this advance, the molecular and atomistic modus operandi of these drugs is yet to come to light. In this study, we employed computational tools to unravel the atomistic interactions and subsequent conformational effects of AMG510 and MRTX849 on the mutant KRASG12C. It was revealed that AMG510 and MRTX849 complexes presented similar total free binding energies, (ΔGbind), of -88.15 ± 5.96 kcal/mol and -88.71 ± 7.70 kcal/mol, respectively. Gly10, Lys16, Thr58, Gly60, Glu62, Glu63, Arg68, Asp69, Met72, His95, Tyr96, Gln99, Arg102 and Val103 interacted prominently with AMG510 and MRTX849. These residues interacted with the pharmacophoric moieties of AMG510 and MRTX849 via hydrogen bonds with decreasing bond lengths at various stages of the simulation. These interactions together with pi-pi stacking, pi-sigma and pi-alkyl interactions induced unfolding of switch I whiles compacting switch II, which could interrupt the binding of effector proteins to these interfaces. These insights present useful atomistic perspectives into the success of AMG510 and MRTX849 which could guide the design of more selective and potent KRAS inhibitors.Communicated by Ramaswamy H. Sarma.

Comparative Dynamic Features of Apo and Bound MDM2 Protein Reveal the Mechanism of Inhibitor Recognition for Anti-Cancer Activity.

BACKGROUND: Mouse Double Minute 2 Homolog (MDM2) oncogenic protein is the principal cellular antagonist of the p53 tumor suppressor gene. Restoration of p53 activity by inhibiting the MDM2-P53 interactions at the molecular level has become the cornerstone of cancer research due to its promising anticancer effects. Natural medicinal products possess various chemical structures and represent an essential source for drug discovery. α-Mangostin (AM) and gambogic acid (G250) are plant-derived compounds that showed inhibitory effects on MDM2-P53 interactions in vitro and in vivo. METHODS: Despite the many clinical studies which performed deeper insight about the molecular understanding of the structural mechanisms exhibited by α-Mangostin and Gambogic acid-binding to MDM2 remains critical. In this study, comparative molecular dynamics simulations were performed for each Apo and bound p53 and MDM2 proteins to shed light on the MDM2-p53 interactions and get a better understanding of the inhibition mechanisms. RESULTS: Results revealed atomistic interaction of AM and G250 within the MDM2-p53 interaction cleft. Both compounds mediate the interaction between the α-helix motifs of the p53 amino-terminal domain, which caused a significant separation between orthogonally opposed residues, specifically Lys8 and Gly47 residues of the p53 and MDM2, respectively. Contrasting changes in magnitudes were observed in per-residue fluctuation on AM and G250 (~0.04 nm and ~2.3 nm, respectively). The Radius of gyration (~0.03 nm and 0.04 nm, respectively), C-alpha deviations (~0.06 nm and 0.1 nm, respectively). The phenolic group of AM was found to establish hydrogen interactions with Glu28 and His96 residues of MDM2. The trioxahexacyclo-ring of G250 also forms hydrogen bond interactions with Lys51 and Leu26 residues of MDM2. CONCLUSION: Utilizing the information provided on the inhibitory binding mode adopted by each compound in this study may further assist in the tailored designs for cancer therapeutics.

Pharmacophore-based Identification of Potential Mutant Isocitrate Dehydrogenases I/2 Inhibitors: An Explorative Avenue in Cancer Drug Design.

BACKGROUND: Heterozygous mutations in the cytoplasmic and mitochondrial isoforms of isocitrate dehydrogenase enzymes 1 and 2 subtypes have been extensively exploited as viable druggable targets, as they decrease the affinity of isocitrate and higher affinity of D-2-hydroxyglutarate, an oncometabolite. OBJECTIVE: Vorasidenib (AG-881) has recently been reported as a promising dual inhibitor of mutant isocitrate dehydrogenase 1 and 2 with the ability to penetrate the blood-brain barrier towards the treatment of low-grade glioma. In order to combat drug resistance and toxicity levels, this compelled us to further investigate this substance as a basis for the creation of potential selective inhibitors of mutant isocitrate dehydrogenases 1 and 2. METHODS: By employing a wide range of computational techniques, binding moieties of AG-881 that contributed towards its selective binding to isocitrate dehydrogenase enzymes 1 and 2 were identified and subsequently used to generate pharmacophore models for the screening of potential inhibitor drugs that were further assessed by their pharmacokinetics and physicochemical properties. RESULTS: AG-881 was identified as the most favorable candidate for isocitrate dehydrogenase enzyme 1, exhibiting a binding free energy of -28.69 kcal/mol. ZINC93978407 was the most favorable candidatefor isocitrate dehydrogenase enzyme 2, displaying a strong binding free energy of -27.10 kcal/mol. ZINC9449923 and ZINC93978407 towards isocitrate dehydrogenase enzyme 1 and 2 showed good protein structural stability with a low radius of gyration values relative to AG-881. CONCLUSION: We investigated that ZINC9449923 of isocitrate dehydrogenase enzyme 1 and ZINC 93978407 of isocitrate dehydrogenase enzyme 2 could serve as promising candidates for the treatment of lower-grade glioma as they cross the blood-brain barrier, and present with lower toxicity levels relative to AG-881.